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We consider quasiparticle propagation in constant-speed-of-sound (iso-tachic) and almost incom- 
pressible (iso-pycnal) hydrodynamic flows, using the technical machinery of general relativity to 
investigate the "effective space-time geometry" that is probed by the quasiparticles. This effective 
geometry, described for the quasiparticles of condensed matter systems by the Painleve-Gullstrand 
metric, generally exhibits curvature (in the sense of Riemann), and many features of quasiparticle 
propagation can be re-phrased in terms of null geodesies. Killing vectors, and Jacobi fields. As par- 
ticular examples of hydrodynamic flow we consider shear flow, a constant-circulation vortex, flow 
past an impenetrable cylinder, and rigid rotation. 

PACS numbers: 02.40.Ky 



I. INTRODUCTION 

The description of many natural phenomena is most 
vividly carried out in terms of hydrodynamics, because 
the concept of a streaming liquid elucidates and helps 
to understand the physical significance and structure of 
an underlying theory Q . In its classical sense 0, 1^ , hy- 
drodynamics describes the motion of a continuum, char- 
acterized by a velocity and density distribution, which 
for a perfect fluid and in the nonrelativistic limit is de- 
scribed by the Euler and continuity equations. It has 
been recognized about twenty years ago by Unruh 0, 
that the propagation of small perturbations on such a 
hydrodynamic background, which is itself governed by 
a continuum version of Newtonian physics, may be cast 
into the form of a "relativistic" scalar wave equation 
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for the velocity potential $ of the perturbations. The dis- 
turbances propagate in an effective space-time with met- 
ric g^v, which is in general curved. The metric g^^y was 
later on shown to be of the Painleve-Gullstrand form , 
originally invented as an alternative to the Schwarzschild 
form of the solution of the Einstein equations for a point 
mass source. With the advent of effective curved space- 
time theories, it became apparent that the Painleve- 
Gullstrand representation of the metric appears in a host 
of such theories. They comprise, besides the conventional 
Euler fluid superfluid ^e-A 0,0, atomic Bose- 

condensed vapors I fll ll?l . and general dielectric (quan- 
tum) matter [Ulllllj 

An interesting and important feature of the Painleve- 
Gullstrand metric is that it continues to give an ap- 
propriate physical description for quasiparticle propaga- 



tion even when the effective space-time possesses a hori- 
zon This occurs because the condensed matter ori- 
gin of the metric in the Painleve-Gullstrand form is the 
spectrum of elementary excitations (quasiparticles) fTHj l. 
which is 'primary. This physical energy spectrum, from 
which the metric is obtained using the fact that for mass- 
less quasiparticles the energy spectrum is 



9 PtJ.Vv = 0, 



(2) 



which must be well-defined and, in particular, real every- 
where in the system. In contrast, for the Schwarzschild 
form of the metric the spectrum reads 
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where rs is the usual Schwarzschild radius and Pr,P± 
are radial and transverse components of the quasipar- 
ticle momentum, respectively. The velocity c plays the 
role of the speed of light and is equal to the sound speed 
for phonons. This "Schwarzschild form" of the spectrum 
exhibits imaginary mode frequencies and consequently 
leads to instability of the condensed matter system if a 
horizon is present, because it has sections of the trans- 
verse momentum p± which result in E'^ < inside the 
horizon. The Painleve-Gullstrand metric, on the other 
hand, gives real frequencies throughout a condensed mat- 
ter system possessing a quasiparticle horizon, which can 
thus be stable. 

The non-equivalence of Schwarzschild and Painleve- 
Gullstrand form of the metric is related to the fact that 
the coordinate transformation relating the Schwarzschild 
solution and the Painleve-Gullstrand representation be- 
comes singular at the horizon 'l4'|. This fact has, inter 
alia, led to the usage of Painleve-Gullstrand co-ordinates 
for investigations of Hawking radiation in the "conven- 
tional" black hole context of gravitational theory [T^It^ . 
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because these co-ordinates are nonsingular through the 
horizon, making the appropriate vacuum definition there 
much simpler. 

The intrinsic characteristics of a curved space-time 
are described in a covariant way by the Riemann tensor 
Our objective in this paper is to describe the 
Riemannian curvature of the effective spaces described 
by the Painleve-GuUstrand metric, in the underlying hy- 
drodynamic terms appropriate to a flowing background 
fluid. We shall focus on two physical situations: quasi- 
particles in flows with a constant speed of sound (iso- 
tachic flows), and quasiparticles in an almost incompress- 
ible (iso-pycnal) hydrodynamic flow. By "almost incom- 
pressible" we mean that we take both the background 
density and the quasiparticle propagation speed relative 
to the medium to be constants, and concentrate on those 
effects that are due to motion of the medium, i.e., its 
velocity distribution. In other words, even if a fluid has 
a constant "refractive index" , focussing and dcfocussing 
effects can be engendered through motion of the fluid. 

As particularly interesting examples we demonstrate 
how the tracks of quasiparticles are distorted by propa- 
gation through a shear flow, a constant-circulation vor- 
tex flow, around an impenetrable cylinder, and how they 
propagate through a rigidly rotating fluid. In a more gen- 
eral context we provide a local definition of "focal length" 
in terms of the Riemann tensor, and show how the afhne 
and "natural" (using the Newtonian background time) 
parameterizations of null geodesies can be related to each 
other. 



II. PAINLEVE-GULLSTRAND CURVATURE IN 
3+1 DIMENSIONS 

In the following discussion the quasiparticle spectrum 
is assumed to be linear in the fluid rest frame for "small" 
quasiparticle momenta, E = c\p\ corresponding to l(2Jl, 
and deviating from linearity for momenta approaching 
the "Planck scale" of the system at hand. In general the 
(3-|-l)-dimensional Painleve-GuUstrand metric |5(| reads 



P r 2 2i 

-- c ~v , gu 
c 
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9u = -- [C - v-\, gu = -- Vi, gij = - Oij. (4) 

That is, the metric has space-time interval 

ds2 = ^ [-cMt^ ^ 5ij{dx' - v'dt){dx^ ~ v^dt)] . (5) 

By special convention, the indices on the 3-velocity are 
always raised and lowered using the flat 3-dimensional 
Cartesian metric so that Vi = u*. 

In the case of irrotational fluid flow (for instance in a 
superfluid outside the cores of the (singular) quantized 
vortices), the d'Alembertian equation (Q) can be derived 
directly from a linearization procedure based on the Eu- 
ler and continuity equations 0] ; the existence and rel- 
evance of the Painleve-GuUstrand effective metric then 
follows as a rigorous theorem. If distributed vorticity 



is present, the situation is more subtle |20j: In hydro- 
dynamics with distributed vorticity one obtains a rather 
complicated system of coupled differential equations, one 
of which contains the d'Alembertian operator (and there- 
fore also contains the effective metric) as a subsidiary 
quantity |20j |. Thus for hydrodynamics with distributed 
vorticity, the effective metric is not the whole story — 
but certainly an important part of the story. In par- 
ticular, if one appeals to the eikonal approximation (in 
this context identical to the WKB approximation) one 
can derive Pierce's approximate wave equation |2ll |. In 
this approximation one can write down the quasiparticle 
spectrum directly in terms of the effective metric ^20] . 

Note that the constant-time hypersurfaces are confor- 
mal to ordinary flat Cartesian space. As long as we are 
interested in quasiparticles that propagate along the null 
cones of this effective metric (that is, quasiparticles mov- 
ing at the speed c relative to the medium), it is permis- 
sible to neglect the overall conformal factor of p/c and 
consider the simplified metric 

gu ^ -[c^ - v^], gu = -Vi, gij = Sij. (6) 

(This is simply the statement that conformal transforma- 
tions leave null curves and, in particular, null geodesies, 
invariant.) The inverse of this simplified metric is 
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Note that the Newtonian time parameter t provides a 
preferred foliation of the spacetime into space-|-time, and 
that this preferred foliation will prove very useful. 

Suppose now that the speed of sound is iso-tachic, in- 
dependent of position and time. Then we can choose co- 
ordinates to set the speed c of linear quasiparticle disper- 
sion equal to unity, a convention adopted in the formulae 
below. The (3-|-l)-dimensional Painleve-GuUstrand met- 
ric then reads 

5tt = -1 + 9ti = -vu 9ij = Sij- (8) 

In general relativistic language the lapse function in the 
ADM formulation is now unity and all the spacetime 
curvature is encoded in the shift function — which here 
describes the physical velocity of the fluid. The inverse 
metric is 



g^- = -i, g- = -v\ g'-' = 5''^ ~ v'^vK (9) 

Turning to the computation of curvature, the 24 inde- 
pendent connection coefhcients read (cf. |2^') 

r*- 



r tt = ViVk Dik = -^{v ■ V)v , 



T'u 



Vj Dij , 



r' jfc — Vi Djk , 

r%t = -dtVi - Vk diVk + ViViVk Dik ■ 



= -dm - \ (S'' - 
= -ViVk Djk + 



(10) 
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Here we have defined the deformation rate and angular 
velocity tensors by 



TrD = divt>, 



5^^u- In a time plus space decomposition 

„ t _ 1 „ t 



= 1, 



ef = , 

j — xJ 

7. ^ l 



(18) 



Thus, for any given vector with components the com- 
ponents in the various frames are related by 



e,/ X^ = {Xf ; Xi) - (Xt + v^X, ; X,) (19) 



(11) and 



The deformation rate is in general relativistic language 
the extrinsic curvature of the constant-time hypersur- 
faces, while the angular velocity tensor is in fluid me- 
chanics language equivalent to the vorticity vector de- 
fined via = e*-''^ ^jk- The above tensors result in the 
unique decomposition of diVj = (V CS) v)ij = Dij + ilij 
into a symmetric and an antisymmetric tensor. 

The components of the Riemann curvature tensor af- 
ford the basic symmetries i?[^,y][pA] — R[pX][t2iy], which 
are supplemented by R[f^i,p\] = and R^[i,p\] = [l^ . 
The Riemann components that need to be calculated 
are thus Rutj, Rijki, and Rujk, the rest follow by the 
(anti-)symmetry properties. A tedious but straightfor- 
ward computation (which follows a variant of the Gauss- 
Codazzi decomposition) yields 



R'ijkl — 
Rtijk = 

Rut] = 



DikDji — DiiDjk , 



(12) 

-diQjk + vi (DkiD^j - DjiD,k) , (13) 
-dtD,, + {Dn + nD\^-{D')^^ 
-VkVk,i] + VkVi {DkiDij - DjkDii) . (14) 

Here we have defined {DVl -f rtD)^j = Dik^lkj + ^ikDkj, 
and similarly [D^^ ^. = DikDkj- 

The appearance and interpretation of the Riemann 
components may be greatly simplified if we consider them 
in an orthonormal, locally Minkowskian tetrad frame 
{e"^}. Greek indices denote the usual spacetime indices, 
Roman letters from the beginning of the alphabet indi- 
cate tetrad indices, while Roman letters from the middle 
of the alphabet denote space indices. Whenever there is 
any chance of confusion, carets on indices are used to in- 
dicate that the components are given in the tetrad frame. 
The tetrad frame {e"^} is defined by 



d^iy — ^ab ^ ^ ^ v • 

In the simplest gauge it is given by 



1 



e t = -V 
The inverse basis satisfies 



51, . 



ab LL h 



(15) 



(16) 



(17) 



Note the use of index placement to distinguish e"^ from 
its inverse Ca^- Hence e°u et^ = '^"b as well as e°^v = 



X" = X^ = (X* ; X') = (X* ; X' - v'X*). (20) 

These index conventions greatly simplify the formulae be- 
low. Calculating the Riemann tensor in the tetrad frame 
gives 



^uDjk , (21) 

(22) 

(D')^^ + {Dn + nD)^^ , (23) 





— DikDji — 


^iifk 


— —difljk , 






Rtitj 





where 



d 



dt+v - V 



(24) 



is the usual convective derivative. The tetrad compo- 
nents Rabcd tell us how a Lagrangian observer moving 
with the fluid perceives the curvature of the effective 
space-time described by the Painleve-GuUstrand metric 
®. 

The components in the tetrad and co-ordinate frames 
are related by 



Ra[3~f5 — ^ a ^ 13 ^ J ^ 5 Rabcd • 

In the tetrad frame, the Ricci tensor 

Rab = R'^acb ~ ^ Rtaib + ^kakb 

has the components 



(25) 
(26) 



Rij 



-^Trr»-Tr(D2), 

^^ikki 

-^(Vxu;), , 

^Riitj + Rkikj 



(27) 



(28) 



dt 



- {Dn + nD)^^ + Tr D, (29) 



where we remind the reader that we have defined the 
vorticity vector 



e^'fe = {rotv)i = (V X v). 



(30) 
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The curvature scalar thus becomes 



III. EXAMPLES 



R 



Rab rj 



ah 



-Rii + Rk 



kk 



2— Tr£) + (TrZ))VTr(D2), (31) 



and contains the trace of the deformation tensor and the 
trace of its square, but not the vorticity. Finally, the 
Einstein tensor takes the form 



Gft 



1 



:R 



i(TrDf -iTr(r>2), 



(32) 
(33) 



- [D,, - SrjTrD) + TrD [D,, - -S,,Tt D 



(34) 



We emphasise that although the Ricci and Einstein 
tensors are nontrivial, and certainly objects of physical 
interest, there is at this level no need for or justification 
for imposing Einstein equations — though these Ricci 
and Einstein tensors are properties of the flow, they are 
not directly related to the stress-energy tensor generating 
that flow and thus the effective space-time curvature ex- 
perienced by the quasiparticles. In superfluids, for exam- 
ple, the "Einstein action" proportional to the curvature 
scalar ()31|l is smaller than the simple kinetic energy of the 
superflow by the factor a^/Z^, where a is the atomic scale 
and I the scale on which the velocity field varies 0], so 
that the "Einstein action" is subdominant in determining 
the velocity field. 

It is sometimes convenient to work with the confor- 
mally invariant, traceless part of curvature. This is given 
by the Weyl tensor . 2 3i] : 

Cabcd = Rabcd + r]a[dRc]b + 'nb[cRd]a + -^R Va[cVd]b, (35) 

where the brackets indicate anti-symmetrization on the 
indices they enclose. This gives 

= Rijki + ^[i Rk]j + + ' (36) 



(37) 



1 d_ 

2dt 



\.. + -(5,„Tr(£)'^ 

u 3 '-^ ^ 



iTr(£>) (^A,- J*y (TrD) 



(38) 



A. General iso-pycnal flows 

Suppose now that the flow is not only iso-tachic (con- 
stant speed of sound) but also iso-pycnal (constant back- 
ground density). This corresponds to an "almost incom- 
pressible" fluid such as water. The major change from 
the previous section is the simplification that comes from 
the continuity equation: 



dt 



V • v = 



Tr D = 0. (39) 



The form of the Riemann tensor is not affected, though 
for the Ricci tensor we now have 



Rti = -Tt{D')., 



Rij — Dij 



{Dfi + nD)^j 

The Ricci scalar simplifies to 

R = Tr(r>2). 



(40) 
(41) 

(42) 
(43) 



Thus the Ricci curvature scalar is positive semidefinite 
for iso-pycnal flows, and vanishes if and only if the defor- 
mation D is zero. 

The Einstein tensor is now 

Git - -iTr(D2), (44) 

Gi, - ^A«,, (45) 

Gij = ^A, - ^%Tr(r»2) _ + UD)^^ . (46) 

Finally the Weyl tensor for iso-pycnal flows reduces to 



^[iRk]j + % + l^^[kSi]p (47) 



^iijk 



ijkl 

difljk + Si[jAvk] , 
1 



1 d 



2 



(AO + OA) . . 



1 



%Tr(A2) . 



(48) 



(49) 



B. Shear flow 

As a first simple example of a nontrivial incompressible 
flow (Tr A — 0), consider the flow with constant shear 



V ^Loo (0,a;,0) 



(50) 



which has both constant deformation D^y = Dy^ = ^luq 
and constant vorticity lOz = coq = 2U,xy = —2U,yx (all 
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other components vanishing) [2J|. The Riemann curva- 
ture components are 



R 



ijki = 4 ^0 - 9u9,k), (51) 



where Oik — ^ki is unity if (ik) — (xy) and zero otherwise. 
The projection operator 



vorticity is non-zero, while in the far field it is the vor- 
ticity that is zero and deformation that is non-zero. The 
Riemann curvature tensor takes the form: 

= -{v\/)D,,-^V.j. (56) 



l^ij — Sij TliTij 



(52) 



More explicitly 



where n = (0, 0, 1) is a unit vector in z direction ensures 
that the curvature has nonzero components only in the 
X and y directions. 

For the Ricci and Einstein tensors 



R-ii = - 0' 



^tyty — ^6 



R 



txty 



{x'-3y'), 
A'j'^xy 



Rii = ^:^^^o^Tt{D^)., 



R 



1 2 
1 2 

1 



Gfi = 0. 



(53) 



Thus the quasiparticles are seen in their effective space- 
time to be moving on a (3-|-l)-dimensional manifold of 
constant scalar curvature, with radius of curvature in- 
versely proportional to the shearing rate t^o- 



Vortex flow of constant circulation 



Riit] 



(57) 



Therefore the Ricci tensor, curvature scalar, and Einstein 
tensor read 



R-ti = 



272 



R = 


27' 


r4 






Gii = 


7^ 
r4 


Gu = 


0. 



72 



(58) 



(59) 



It is mildly amusing to note that the vortex geometry 
is uniquely determined by the cylindrical symmetry plus 
the equation Gah 8ah (not r\ah)- 



A somewhat more interesting case is the constant- 
circulation flow in the x-y plane 



D. Streaming motion past a cylinder 



7a; 



72/ 



x-^ -I- 



x^ -\- 



(54) 



appropriate to a vortex flow well outside the central core, 
where the circulation is f v ■ ds = 2t:^. In this case you 
would not want to trust the geometry for r < Tc = 7 
because at r — Tc the flow goes supersonic. This flow has 



n,. 



A - ^yy 



0. 



= A, 



(55) 



The most complex flow we discuss here is provided 
by the two-dimensional streaming motion from right to 
left past a cylinder of radius a. According to the circle 
theorem ^] , the complex velocity potential of such a flow 
is given by 



w = U\Z+ — 



(60) 



where Z = x + iy and U is the velocity at infinity in 
negative x direction. This results in the flow 



-C/ ( 1 -t- a2 ^^-j^ V Vy = 2Uxy^. (61) 



Note the "duality" between the vortex core and the far 
field. In the core the deformation rate is zero and the 



The velocity at infinity is restricted to be (7 < 1/2, for 
the maximal velocity on the cylinder surface to be less 
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than the speed of sound. The formulae for deformation 
and vorticity (which is identically zero for this flow) read 





2Ua^ , 




Dxy 




= ^6 y( 


A. 






= 0. 



D 



yx 



(62) 

The Riemann components show that the flow past a 
cylinder, due to its reduced symmetry, yields a more com- 
plicated space-time geometry for quasiparticles than the 
vortex flow: 

1, 



Riyiy = det£) = --Tr(D^) 



= -{v- V)D,j + Vij det D , (63) 
where the last line reads more explicitly 



[a^(y2-5:r2)-|-3(x4-6:rV + /)], 



R. 



2„2 



tyty 



^ixiy 



-X 



(64) 



These latter components show that the "circulation" Ua? 
is not the only relevant parameter of the flow, in contrast 
to the constant-circulation vortex case, as we may expect 
from the reduced symmetry of the flow past the cylinder. 
The curvature scalar 



R 



(65) 



decays much more quickly with distance from the cylin- 
drical object than the curvature of the vortex flow, Eq. 

(EH- 



E. Rigid rotation 

The simplest example of a nontrivial incompressible 
flow (Tr£> — 0) is pure rotation v — n{~y,x,0), which 
has zero deformation Dij = 0, and constant vorticity 
LUz = loq = 2Qxy — ~2Qyx — 2Q (all other components 
vanishing) . This flow is appropriate for instance deep in- 
side the core of a vortex where the fluid effectively rotates 
as a "rigid" body. (In ordinary fluids this happens be- 
cause viscosity dominates in the core; in superfluids there 
is a more dramatic effect in that the superfluid goes nor- 
mal close enough to the core.) Also note that the core 
has a maximum size given by = 1, that is Vc = S/wq. 



For the rigid rotation flow it is easy to see that the 
Riemann curvature tensor is identically zero, either (1) 
by brute force application of the above formulae, or more 
subtly (2) by going to a rotating frame (of angular ve- 
locity Q — 2uJo) in which the velocity is identically zero, 
evaluating the Riemann tensor there (where it is bla- 
tantly zero), and transforming back to the rotating frame. 
Although the Riemann tensor is identically zero, there is 
interesting physics going on: The fact that pure rotation 
leads to zero Riemann curvature is ultimately responsible 
for the fact that equations (|12|l and (|21|l do not contain 
any terms quadratic in f2, a result that otherwise has to 
be simply asserted based on explicit calculation. 

Additionally, we emphasise that even though the Rie- 
mann tensor is zero, the Christoffel symbols are definitely 
not zero. Indeed 



1 



(66) 
(67) 



These two portions of the Christoffel symbols are of 
course simply representing the centrifugal and Coriolis 
pseudo- forces. All other components are zero. 

A further (approximate) example of such a flow is en- 
countered if one considers the coarse-grained flow in- 
duced by a lattice of vortices (inflnite) lattice 
rotates as if it were a solid body, with a vortex density 
— ^/n^ prescribed by the rotation velocity f2 and 
the circulation 27r7, assumed to be equal for each indi- 
vidual vortex. For the vortex lattice, it follows from the 
vanishing of the Riemann curvature that a collimated 
quasiparticle beam can pass a (sufficiently dilute) lattice 
without (on average) being deflected. 



IV. GEODESIC DEVIATION 

An invariant measure of the strength of a flow pattern 
as regards its influence on quasiparticle motion may be 
defined to be the value of the curvature scalar R (x s~'^ 
at a certain given distance s from the fiow-generating 
object (cf. Fig. ^ illustrating the generic situation of 
flow past an object placed in a homogeneous stream). 
Among the flows discussed in the previous section the 
shear flow is strongest in that sense (because the "flow 
generating object" is covering all space, k = 0), followed 
by the vortex flow {k ~ 4) and the flow past the cylinder 
(k = 6). Finally rigid rotation, which has zero R and is 
"flat" (k = oo). It is the simplest conceivable nontrivial 
(i.e., inhomogeneous) flow with the property of having 
all Rabcd equal to zero. 

A nonvanishing Riemann tensor leads to tidal (rela- 
tive) acceleration of nearby geodesies, described by the 
Jacobi equation of geodesic deviation for quasiparticles: 



~dA2~ 



+ R'^fSjsu'^n^u^ 



0. 



(68) 
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FIG. 1: The quasiparticle geodesic deviation at a distance 
vector s caused by an object placed in a flow with velocity Voo 
at infinity (the generic case of the situation in section Hit Dl l. 



The above relation gives the covariant relative accelera- 
tion of two nearby geodesies, with null tangent vectors 
u separated by the displacement vector n, and with the 
geodesies afSnely parametrized by A. (At this stage all we 
need to know is that use of an "affine parameter" simpli- 
fies many formulae; in the following section we will derive 
a relationship between the afhne parameter and physical 
Newtonian time t.) The fact that the constant time slices 
of the metric Q are conformally identical to flat Carte- 
sian space in three dimensions, entails that the space- 
time curvature of the quasiparticle world is reflected in a 
relative acceleration of quasiparticle rays in the Newto- 
nian lab world of nonrelativistic hydrodynamic flow. 

Consider a family of geodesies in the z-direction, with 
tangent vector = (u*,ii*, 0,0) and a purely spacehkc 
separation in the z/-direction n — (0, 0, Jy, 0). We then 
have 



dA2 



{( 



R 



ytyt 



Ryxyx) (u ) I 



0. 



(69) 



This can be viewed as a parametrically driven harmonic 
oscillator (driven in the affine parameter A), with "fre- 
quency" 



n{X) 



ytyt 



Ryxyx • 



(70) 



Physically this means that by looking at the components 
of the Riemann tensor we can see if the effective geome- 
try locally acts as a focussing lens [corresponding to r2(A) 
real] or as a diverging lens [corresponding to Q{X) imag- 
inary]. Since (in the focussing case, and assuming a rea- 
sonably uniform medium) two initially parallel geodesies 
will focus down to a point after an elapse of affine param- 
eter 5A = 7r/ri(A), the corresponding local focal length is 
(in physical distance units) given by 



/ 



local 



±7r 



-^iyiy + ^xyxy \ 



(71) 



Note the strengths and weaknesses of this concept — 
it provides a local position and orientation dependent 
notion of focal length appropriate for nearly parallel 



geodesies (nearly parallel quasiparticles; so one is au- 
tomatically working "on axis" and ignoring "spherical 
abberation"), but this definition of f^°'^'^^ does in gen- 
eral not provide significant global information. If the 
Riemann tensor is strongly inhomogencous, varying on 
length scales significantly smaller than f^°'^^^^ then this 
concept of local focal length is not particularly useful. 
In particular, in the vortex geometry of ,26] , with flow 
the focussing effect we had in mind was a global 
effect due to quasiparticles passing by opposite sides of 
the vortex core, with impact parameter b — this is not 
a situation that can be described by the Jacobi equa- 
tion. The global result obtained there for / — /gi°t>ai _ 
{2b^/3Trrl)[l + 0{rc/b)], is not the local /'"'^^i defined 
above. Indeed two initially parallel quasiparticles pass- 
ing by on the same side of the vortex core will be driven 
apart from each other by geodesic deviation — it is this 
effect that leads to the "cylindrical abberation" of the 
lens discussed in |26j . 

A case where the local focal length does acquire global 
meaning is the shear flow (|50|l . for which the focal length 
(|7T|l becomes a constant: 



ihear 



(72) 



The focal length is in this case bounded by the atomic 
length scale itself, simply due to the requirement that 
the concept of hydrodynamics makes sense. This fur- 
ther strengthens the notion of the shear flow being the 
strongest possible flow as regards its influence on quasi- 
particle motion, because any other flow has more strin- 
gent bounds on the global f. 

One useful refinement of the local focal length concept 
introduced in equation H71|l is to consider null geodesies 
(quasiparticle paths) propagating in an arbitrary unit di- 
rection u and then use indices M and N to denote the 
two spatial directions perpendicular to H. Then the local 
focal length can be generalized to a 2 x 2 matrix 



^local 
JMN 



±7r 



\/\\RiMiN + RiMjN Ul\\ 



(73) 



The square root and inverse is to be taken in the ma- 
trix sense, and the two eigenvalues of /a/ at are the two 
principal focal lengths along the direction u. If these 
eigenvalues differ it is a signal of astigmatism. 



V. NON-AFFINE PARAMETERIZATION OF 
NULL GEODESICS 

While the use of affine parameters for null geodesies is 
standard in general relativity, it should be borne in mind 
that in the present Painleve-GuUstrand context there is a 
preferred temporal foliation provided by the Newtonian 
time parameter t. It is worth the technical bother of 
using the non-affine parameterization in terms of t here 
in order to make aspects of the physics clearer. 



8 



In general, we know that along any null geodesic there 
will be some relationship between afhne parameter A and 
Newtonian time t. For instance we can assert 



dA = exp[C(t)] dt. 



(74) 



In the afhne parameterization the geodesic equation for 
a null curve is just 



0; u'' = — 



dX' 



If we choose a non-affine parameterization 

dr" 

u^V^r = at) w"; = 
The geodesic equation becomes 



d'^x^' 
~dt2~ 



a/3 



dx" dx^ 

"dT ~dr 



da;'' 
"dT 



(75) 



In this form it is clear that the physical acceleration of 
the quasiparticle is related to gradients in the Painleve- 
GuUstrand metric. It is extremely useful to derive an 
explicit relationship between the affine parameter A and 
the physical Newtonian time t. To do this let's start 
with the notion of a stationary geometry (technically: 
there exists a timelike Killing vector; colloquially: a time- 
independent geometry) . The timelike Killing vector takes 
the form 

K^ = {1;6); K^ = {-[l-v^]; -v). (76) 
The tangent vector to the null geodesic is denoted 

da;' 



da;'^ 
"dA" 



dt 
dA 



1; 



dt 



(77) 



It is a standard theorem that the 3+1 inner product be- 
tween a geodesic tangent vector and a Killing vector is 
conserved, as long as the geodesic is affinely parameter- 
ized. Thus 



2 da; 
dt 



~ constant. (78) 



On the other hand, because u^^ is a null vector 



1 



2v- 



dx 
'dt 



dx 



dt 



= 0. 



(79) 



Eliminating between these two equations, we can normal- 
ize in such a way that 



dt r 

— = exp[-C(t)] = 



1 _ + 



da; 
dt 



That is 



at) = In 



da; 
lit 



(80) 



(81) 



If the fluid is not moving, then v = and |da;/dt| = 
1 so t (X A. If the fluid is moving we simply have to 
live with this position-dependent factor relating the affine 
parameter A (in terms of which the geodesic equations 
are most easily written down) to the Newtonian time 
parameter t (in terms of which the physical acceleration 
is most easily calculated). 

In a similar manner, the Jacobi equation can be rewrit- 
ten as 



dt2 



Dt 



dt 



0. 



(82) 



While this looks somewhat messier than the affinely pa- 
rameterized Jacobi equation Ht)8|) , the physics is the same. 
In particular if we start with two initially parallel null 
geodesies [Dn/dt = at t = 0), and assume a locally 
homogeneous medium, we are led to the same notion of 
local focal length as discussed in the previous section. 



VI. DISCUSSION 

We have shown how the generation of curved Ricman- 
nian space-time geometries for quasiparticles is possible 
based purely on the velocity pattern of a nonrelativistic 
flow. Conversely, one might conceive of solving for a flow 
fleld from a given space-time geometry. This is a highly 
nonlinear problem, as becomes obvious from the relations 
(|21|I - H23|) . It is, however, certainly no more nonlinear or 
complicated than solving the Einstein equations of gen- 
eral relativity themselves. While the Painleve-GuUstrand 
geometry discussed here does not provide us with the 
most generic case (remember that the constant time sur- 
faces are (conformally) flat; for generalizations allowing 
for more general space-time metrics see [lOj ) . it shows 
that the underlying kinematical structure of a curved 
space-time can in principle be perfectly nonrelativistic. 
The dynamical identification of this effective geometry 
with general relativity, i.e., imposing the Einstein equa- 
tions, is a more advanced step f2|, but is possible in prin- 
ciple as well. 

There are several generalizations of the current anal- 
ysis that would be of interest: (1) If the quasiparticle 
propagation speed (c, local speed with respect to the 
background medium) is varying then the geometry ex- 
hibits "index gradient" effects in addition to effects gen- 
erated by the motion of the medium. While technically 
straightforward, the relevant calculations of the Riemann 
tensor are computationally messy and the physical inter- 
pretation is not so clear (unless the medium is completely 
at rest; in which case one recovers standard "index gra- 
dient" physics). (2) If the density varies from place to 
place, then it is necessary to distinguish the "geomet- 
rical quasiparticle" regime (the analogue of geometrical 
optics) from the "wave quasiparticle regime" (the ana- 
logue of wave optics). In the geometrical approximation 
the results of the present paper can be carried over; in 
the wave regime one needs to carry out an analysis in 
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terms of Green functions and wave equations; the entire 
armoury of quasiparticle trajectories as null geodesies of 
the effective metric breaks down and must be replaced 
by a more fundamental wave description. 

In summary: The use of pseudo-Riemannian geome- 
try has important applications well beyond the confines 
of general relativity. In particular quasiparticle propa- 
gation in condensed matter systems can often be char- 
acterized in terms of an "effective" spacetime geometry; 
most easily described in Painleve-GuUstrand form. If the 
background medium is a fluid, then the Ricmann curva- 
ture (and Christoffel symbols, etc.) can be calculated in 



terms of shear (deformation) and vorticity of the fluid. 
Ultimately this analysis relates the focussing and deflec- 
tion of quasiparticles to the properties of the fluid flow. 
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